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PI:  Nan  Guang  Chen 


Introduction: 

Breast  cancer  is  the  number  two  killer  of  women  in  United  States.  Mammography  and 
ultrasound  are  the  most  common  imaging  modalities  for  breast  cancer  screening  and 
diagnosis.  However,  X-rays  see  very  low  contrasts  for  human  soft  tissues.  It’s  hard  to 
differentiate  tumors  from  other  lesions  such  as  cysts  on  a  mammogram.  Ultrasound  can 
detect  and  resolve  breast  lesions  a  few  millimeters  in  size,  and  can  distinguish  cysts  from 
solid  lesions.  Nonetheless,  the  overlapping  acoustic  characteristics  of  benign  and 
malignant  tumors  result  in  a  low  specificity  of  ultrasound.  Currently,  only  20%  patients 
sent  to  biopsy  have  malignant  tumors,  and  80%  patients  suffer  from  unnecessary  biopsy. 
The  purpose  of  this  research  is  to  develop  a  combined  ultrasound  and  near-infrared 
diffusive  imager  that  can  significantly  improve  the  specificity  for  breast  cancer  detection 
and  reduce  the  unnecessary  biopsy  rate.  Functional  imaging  with  near-infrared  light  is 
based  on  diffusive  propagation  of  photons  inside  turbid  media  such  as  breast  tissues.  It 
has  been  reported  that  this  new  imaging  modality  has  the  potential  to  differentiate  breast 
tumors  from  normal  surrounding  tissues  through  determination  of  optical  properties  at 
different  wavelengths.  Some  physiological  parameters,  such  as  blood  volume, 
hemoglobin  oxygenation  level,  and  water  concentrations  can  be  retrieved  from  optical 
measurements.  However,  the  diffusive  optical  imaging  method  alone  cannot  readily  be 
applied  to  clinical  diagnosis  because  of  its  low  spatial  resolution  and  some  difficulties  in 
image  reconstruction.  Combination  of  ultrasound  and  near-infrared  diffusive  imaging 
benefits  from  high  spatial  resolution  of  ultrasound  and  functional  contrasts  from  optical 
imaging.  In  addition,  a  priori  morphological  information  about  embedded  lesions  from 
ultrasound  can  help  improve  the  accuracy  of  optical  image  reconstruction.  In  this 
research,  we  have  investigated  the  probe  design  involving  arrangement  of  ultrasound 
transducer  and  optical  light  guides,  and  optimization  of  source  and  detector  numbers.  A 
frequency-domain  and  a  DC  diffusive  optical  imaging  systems  have  been  built  and 
calibrated.  Phantom  studies  have  been  carried  out,  and  new  image  reconstmction 
algorithms  have  been  proposed  and  validated.  Our  research  results  have  been  reported  in 
a  series  of  journal  papers  and  conference  presentations.  In  the  next  step,  we  will  use  our 
combined  imager  in  Hartford  Hospital  and  UConn  Health  Center  and  gather  clinical 
experimental  results.  We  will  look  into  further  problems  emerging  from  clinical 
situations  and  improve  our  imaging  systems  as  well  as  image  reconstmction  algorithms. 
Correlations  between  combined  ultrasound/optieal  images  with  biopsy  results  will  be 
established  after  clinical  studies. 

Body: 

The  optical  sources  and  detector  numbers  and  their  deployment  on  a  probe  are  among 
critical  issues  concerning  diffusive  optieal  imaging.  Because  of  the  diffusive  nature  of 
near-infrared  light  migrating  through  human  soft  tissues,  more  light  sources  and  sensors 
do  not  necessarily  provide  more  independent  information.  On  the  other  hand,  the  system 
eost  and  data  acquisition  time  increase  with  the  numbers  of  sources  and  detectors.  We 
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conducted  computer  simulation  to  find  the  optimal  number  and  the  best  arrangement.  In 
our  simulation,  a  total  of  364  optical  source-detector  pairs  were  deployed  uniformly  over 
a  9  cm  X  9  cm  probe  area  initially.  Then  the  total  pairs  were  reduced  gradually  to  60  in 
experimental  and  simulation  studies.  For  each  source-detector  configuration,  three- 
dimensional  images  of  a  1-cm-diameter  absorber  of  different  contrasts  were  reconstmcted 
from  the  measurements  made  with  a  frequency-domain  system.  The  results  have  shown 
that  more  than  160  source-detector  pairs  are  needed  to  reconstruct  the  absorption 
coefficient  to  within  60%  of  the  true  value  and  appropriate  spatial  and  contrast  resolution. 
However,  the  error  in  target  depth  was  more  than  1  cm  in  all  source-detector 
configurations.  With  the  a  priori  target  depth  information  provided  by  ultrasound,  the 
accuracy  of  the  reconstructed  absorption  coefficient  was  improved  by  15%  and  30%  on 
average.  The  speed  of  reconstruction  was  improved  by  ten  times  on  average.  As  the 
ultrasound  transducer  needs  to  be  placed  in  the  middle  of  a  probe,  we  have 
experimentally  evaluated  the  effects  of  missing  optical  sensors  in  the  middle  of  the 
combined  probe  upon  the  accuracy  of  the  reconstructed  optical  absorption  coefficient, 
and  assessed  the  improvements  of  reconstructed  absorption  coefficient  with  the  guidance 
of  the  co-registered  ultrasound.  The  results  have  shown  that  when  the  central  ultrasound 
array  area  is  in  the  neighborhood  of  2  x  2  cm^,  which  corresponds  to  the  size  of  most 
commercial  ultrasound  transducers,  the  optical  image  quality  is  not  affected.  The  results 
also  shown  that  the  iterative  inversion  algorithm  converges  very  fast  with  the  guidance  of 
a  priori  target  spatial  distributions. 

We  have  constructed  a  near-infrared  frequency-domain  imaging  system.  This  system 
features  multiple  transmission,  parallel  reception,  and  dual  wavelengths.  On  the 
transmission  part,  12  pairs  of  dual  wavelength  (780nm  and  830nm)  laser  diodes  are  used 
as  light  sources,  and  their  outputs  are  amplitude  modulated  at  140.000  MHz.  To  reduce 
noise  and  interference,  an  individual  driving  circuit  is  built  for  each  diode.  As  laser  diode 
works  in  series,  a  control  board  that  interprets  instructions  from  a  PC  is  used  to 
coordinate  operations  of  associated  components.  When  one  transmission  channel  is 
selected,  it  turns  on  corresponding  driving  circuit  so  that  a  DC  driving  current  can  be  set 
up  for  the  diode.  At  the  same  time,  a  selected  signal  is  sent  to  a  RF  switching  unit,  which 
distributes  RF  signal  to  the  right  channel  to  modulate  the  optical  output.  On  the  reception 
part,  8  Photo  Multiplier  Tubes  (PMTs)  are  employed  to  detect  diffusely  reflected  light 
from  turbid  media.  Each  PMT  is  housed  in  a  sealed  aluminum  box,  shielding  both 
environment  lights  and  electromagnetic  fields,  while  an  optical  fiber  (3  mm  in  diameter) 
couples  NER  light  from  the  detection  point  to  the  reception  window  of  the  PMT.  The 
electrical  signal  converted  from  optical  input  is  generally  very  weak  and  rather  high  in 
frequency,  so  high  gain  amplification  and  frequency  transform  are  necessary  before  it  can 
be  sampled  by  an  A/D  board  inside  the  PC.  We  have  built  8  parallel  heterodyne 
amplification  channels  so  as  to  measure  the  response  of  all  detectors  simultaneously,  thus 
to  reduce  the  data  acquisition  time.  Each  amplification  channel  consists  of  a  RF 
amplifier  (40dB),  a  mixer  where  RF  signal  (140.000  MHz)  is  mixed  with  a  local 
oscillator  (140.020  MHz)  and,  a  band  pass  filter  centered  at  20  KHz,  and  a  low  frequency 
amplifier  of  30dB.  The  heterodyned  two-stage  amplification  scheme  helps  suppress  wide 
band  noises  efficiently.  Eight  detection  signals  are  sampled,  converted,  and  acquired  into 
the  PC  simultaneously,  in  which  Hilbert  transform  is  used  to  compute  the  amplitude  and 


5 


phase  of  each  channel.  The  entire  data  acquisition  is  less  than  15  seconds,  which  is  fast 
enough  to  acquire  data  from  patients. 

A  DC  system  has  also  been  built  for  comparing  its  performance  with  that  of  a  frequency- 
domain  system.  The  DC  system  has  6  pairs  of  dual  wavelength  sources  (780  nm  and  830 
nm).  There  are  12  detection  light  guides  that  collect  photons  from  12  positions  and  are 
coupled  to  a  PMT  through  a  mechanical  multiplexer.  The  multiplexer  is  implemented  by 
a  rotating  plate  driven  by  a  stepper  motor.  Although  the  total  data  acquisition  speed  is 
slow  (about  1  minute  for  a  complete  scan),  the  DC  system  is  far  more  compact  than  the 
frequency  domain  system  and  costs  much  less.  We  have  tested  and  calibrated  both 
systems  with  phantoms. 

To  avoid  the  problems  related  to  iterative  image  reconstruction  algorithms,  we  proposed 
a  novel  non-iterative  method  using  minimal  a  priori  co-registered  ultrasound 
information.  Small  absorbing  targets  embedded  in  a  homogeneous  background  are 
described  approximately  in  terms  of  their  monopole,  dipole,  and  quadrupole  moments. 
With  an  approximate  estimation  of  center  locations  of  these  absorbers  from  ultrasound 
images,  we  show  in  simulations  that  the  reconstruction  accuracy  in  absorption  coefficient 
exceeds  80%  if  the  noise  level  is  less  than  0.2%.  We  also  demonstrate  experimentally 
that  the  accuracy  can  be  improved  by  using  additional  ultrasound  volume  information 
even  for  a  noise  level  as  high  as  1.5%.  In  the  moment-based  method,  we  assume  that  the 
lesions  are  isolated  and  are  imbedded  in  a  homogeneous  background.  This  assumption  is 
quite  true  for  more  homogeneous  fatty  breasts  and  may  not  hold  for  dense  breasts  that 
consist  of  both  glandular  tissue  and  fat.  In  the  latter  case,  we  could  segment  ultrasound 
images,  identify  tissue  types  and  estimate  background  optical  properties  of  different 
tissues  in  the  reconstruction  as  well.  Both  acoustic  and  optical  contrasts  exist  in  tumors 
but  the  sensitivities  of  these  two  modalities  may  be  different.  Therefore,  correlation 
between  acoustic  and  optical  heterogeneities  remains  unknown  at  the  current  stage  and 
will  be  answered  by  our  future  clinical  studies. 

Key  research  accomplishments. 

•  Built  a  frequency-domain  near-infrared  diffusive  imaging  system,  which  has  12 
pairs  of  dual-wavelength  sources  and  8  detectors. 

•  Built  a  DC  near-infrared  diffusive  imaging  system,  which  has  6  pairs  of  dual¬ 
wavelength  sources  and  12  detector  positions  coupled  to  1  PMT  by  a  mechanical 
multiplexer. 

•  Tested  and  calibrated  both  the  frequency-domain  and  DC  systems  with  phantoms. 

•  Found  the  optimal  numbers  of  sources  and  detectors  from  computer  simulations. 

•  Conducted  phantom  studies  with  combined  ultrasound/optical  imager. 
Investigated  the  effect  of  central  ultrasound  transducer  on  the  optical  image 
quality.  Demonstrated  that  a  priori  target  information  from  ultrasound 
significantly  improved  the  accuracy  of  reconstructed  values  of  optical  properties. 

•  Propose  a  novel  diffusive  image  reconstruction  method  based  on  lower  order 
moments  of  embedded  targets  inside  homogeneous  turbid  media. 
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Reportable  outcomes: 

Journal  papers: 

[1]  N.  G.  Chen  and  Q.  Zhu,  “Characterization  of  small  absorbers  inside  turbid  media,” 
Optics  Letters  ,  Vol  27,  No.  4,  pp  252-254  (2002). 

[2]  Q.  Zhu,  N.  G.  Chen,  ,  P.  Y.  Guo,  ,  S.  K.  Yan,  and  D.  Q.  Piao,  ,  “Near  infrared 
diffusive  light  imaging  with  ultrasound  localization,”  OSA  Optics  and  Photonics  News, 
Optics  in  2001,  12:31,  December  2001 

[3]  N.  G.  Chen,  P.  Guo,  S.  Yan,  D.  Piao,  and  Q.  Zhu,  “Simultaneous  near  infrared 
diffusive  light  and  ultrasound  imaging,”  Applied  Optics,  vol.  40,  No.  34,  pp  6367-6280, 
2001. 

[4]  Q.  Zhu,  N.  G.  Chen,  X.  H.  Ding,  D.  Piao,  and  P.  Guo,  “Design  of  near  infrared 
imaging  probe  with  the  assistance  of  ultrasound  localization,”  Applied  Optics,  vol.  40, 
No.  19,  pp  3288-3303,2001. 

Presentations: 

[1]  N.  G.  Chen,  and  Q.  Zhu,  “Nonlinear  correction  method  for  characterizing  small 
absorbers  in  turbid  media,”  OSA  Biomedical  Topical  Meeting,  Miami  Beach,  March 
2002. 

[2]  Minming  Huang,  Tuqiang  Xie,  Nan  Guang  Chen,  and  Quing  Zhu,  “NIR  imaging 
reconstruction  with  ultrasound  guidance:  Finite  element  method,”  OSA  Biomedical 
Topical  Meeting,  Miami  Beach,  March  2002. 

Funding  applied  for: 

[1]  Donaghue  Medical  Research  Foundation:  Co-PI,  “3D  simultaneous  Ultrasound  and 
Near  Infrared  Diffusive  Imaging  for  Breast  Cancer  Detection;” 

[2]  NIH:  Co-PI,  “Near  Infrared  Diffused  Light  Imaging  with  Ultrasound  Guidance.” 

Conclusions: 

I  have  successfully  finished  the  Task  1  specified  in  the  statement  of  work  of  my  proposal. 
A  DC  and  a  frequency-domain  near-infrared  diffusive  imaging  systems  have  been  built, 
tested,  and  calibrated.  Phantom  studies  have  demonstrated  the  advantages  of  the 
combined  ultrasound/optical  imager.  A  novel  image  reconstruction  algorithm  has  been 
proposed  and  implemented.  I  am  preparing  for  the  clinical  study  that  will  eventually  test 
the  effectiveness  of  this  new  technique  for  breast  cancer  detection. 
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We  propose  a  novel  noniterative  near-infirared  diffusive  image  reconstruction  method  that  uses  minimal 
a  priori  co-registered  ultrasound  information.  Small  absorbing  targets  embedded  in  a  homogeneous 
background  are  described  approximately  in  terms  of  their  monopole,  dipole,  and  quadrupole  moments.  With 
an  approximate  estimation  of  the  center  locations  of  these  absorbers  Scorn  ultrasound  images,  we  show  in 
simulations  that  the  reconstruction  accuracy  of  the  absorption  coefficient  exceeds  80%  if  the  noise  level  is 
less  than  0.2%.  We  also  demonstrate  experimentally  that  the  accuracy  can  be  improved  by  use  of  additional 
ultrasound  volume  inforniation  even  for  a  noise  level  as  hi^  as  1.5%.  G  2002  Optical  Society  of  America 
OCIS  codes:  100.3190,  170.3010,  170.3880,  290.1990. 


Functional  imaging  with  near-infrared  (NIR)  light 
has  found  potential  applications  in  many  areas, 
such  as  breast  and  brain  lesion  detection  and  diag¬ 
nosis.^  Recently,  a  combination  of  NIR  imaging 
with  other  imaging  modalities,  such  as  ultrasound  or 
magnetic  resonance  imaging,  has  shown  promising 
results*'^  in  providing  complementary  contrasts  and 
overcoming  NIR  reconstruction  problems  related  to 
intensive  light  scattering. 

In  this  Letter  we  introduce  a  novel  reconstruction 
algorithm  for  NIR  diffusive  imaging  that  uses  approxi¬ 
mate  target  center  locations  estimated  from  co- 
r^fistered  ultrasoimd.  Oiu*  new  method  is  based 
on  estimation  of  major  characteristics  of  isolated 
small  absorbers.  These  characteristics  are  monopole, 
dipole,  and  quadrupole  moments.  Higher-order  mo¬ 
ments  have  negligible  effects  in  characterizing  the 
absorbers  and  are  ignored.  It  is  shown  in  this  Letter 
that  measurements  of  diffusive  photon  density  waves 
cannot  readily  achieve  a  signal-to-noise  ratio  that  is 
good  enough  for  reconstructing  the  detailed  shape 
information  of  targets.  However,  the  quadrupole 
moment  can  provide  an  approximate  extension  of  the 
target,  which  is  necessary  for  distributing  the  integral 
absorption  (the  monopole)  to  an  appropriate  target  re¬ 
gion.  In  situations  in  which  the  measurement  system 
suffers  from  unexpected  noise,  accurate  estimation  of 
the  quadrupole  moment  might  not  be  possible.  The 
target  volume  estimated  from  co-registered  ultrasound 
images  can  be  used  for  distributing  the  integ^ral  ab¬ 
sorption  inside  the  target  volume.  We  have  conducted 
a  series  of  simulations  at  different  noise  levels  to 
evaluate  the  performance  of  this  new  method.  We 
have  also  conducted  experiments  using  our  combined 
NIR-ultrasound  imager®  to  test  the  algorithm. 

If  it  is  assumed  that  there  are  only  a  few  small  inho¬ 
mogeneous  taigets  embedded  in  a  homogeneous  back¬ 
ground,  the  distribution  of  the  absorption  coefficient 
can  be  expressed  as 

flair)  =  fla^  +  A/ia(r) .  (1) 

The  background  absorption  coefficient,  fia^,  can  be  as 
small  as  0.02  cm“‘  for  normal  breast  tissues,  whereas 
Afia  can  be  well  beyond  0.1  cm~*  for  tumors  because 
of  the  blood.  In  the  Bom  approximation,  the  photon 
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density  wave  is  decomposed  into  an  incident  wave  and 
a  scattering  wave: 

r)  =  4>i„c(6>,  r)  +  r) .  (2) 

The  incident  wave  is  the  solution  of  a  homogeneous  dif¬ 
fusion  equation,  and  the  scattering  wave  can  be  given 
approximately  by  the  linear  perturbation  theory: 

^^^^iacio>y)G(r,r')d^r'.  (3) 

In  this  equation,  Gir,rO  is  the  Green’s  function, 
and  D  =  1/3/i,'  is  inversely  proportional  to  the  re¬ 
duced  scattering  coefficient,  /*,'.  In  the  context  of 
early  stage  Ibreast  cancer  detection,  we  can  further 
assume  that  those  heterogeneities  are  confined  to  a 
few  insolated  regions.  Then,  the  scattering  wave  can 
be  expanded  aroimd  each  of  N  target  centers  as 

H  Z  [M-'WCro”)  -I-  D"  •  VWfro”) 

-I-  Q”  •  VVW(ro*')/2  -f-  0(a®)] .  (4) 

Here  ro"  is  the  center  of  the  yth  target,  and  Wfro”)  = 
4*inc(<«».ro’')G(r,ro’’)  is  the  weight  fimction.  For  the 
I'th  tai^et,  M",  D',  and  Q’’  are  its  monopole,  dipole,  and 
quadrupole  moments,  respectively.  Terms  beyond  the 
second  order  of  target  dimension  a  have  negligible  ef¬ 
fects  in  characterizing  the  absorbers  and  are  neglected 
in  our  model.  For  the  I'th  target,  its  monopole,  dipole, 
and  quadrupole  moments  are  scalar,  vector,  and  sec¬ 
ond-order  tensor,  respectively,  and  they  are  given  as 

M"  =  f  Afiair')d^r' ,  D"  =  f  A/ta(r')r'd®r', 

Jv,  Jv, 

0"=  f  A/ia(r')r'r'dV.  (5) 

Jv, 

The  integrations  are  over  a  small  isolated  region  V,. 

A  multiple-source,  multiple-detector  configuration 
is  typical  for  frequency-domain  diffusive  imaging 
systems.  We  are  using  12  sources  (780  nm)  that  are 
amplitude  modulated  at  140  MHz  and  eight  detectors 
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in  our  simulations  as  well  as  in  the  experiments. 
These  sources  and  detectors  are  deployed  pn  an  ab¬ 
sorbing  plane,  which  simplifies  boundary  conditions. 
The  reflection  mode  with  a  semi-infinite  geometry 
is  used  for  both  simulation  and  experiment.  The 
measured  scattering  wave  is  related  to  the  target 
moments  by  the  following  equation: 

[^i--4mf  ...  Q^f.  (6) 

An  element  of  11  is  either  the  weight  function  for  one 
target  or  its  derivatives  up  to  the  second  order. 

The  inverse  problem  is  to  retrieve  the  character¬ 
istics  of  those  embedded  absorbers  from  measured 
photon  density  waves  on  the  surface.  Based  on  our 
simpliHed  forward  model,  it  is  necessary  to  have 
certain  information  about  absorbers,  such  as  the  num¬ 
ber  of  targets  and  their  center  positions,  in  advance. 
The  ultrasound  system  of  our  combined  imager  can 
provide  sudi  information.  Once  the  initial  center 
locations  have  been  specified,  the  weight  function, 
together  with  its  first  and  second  derivatives,  can  be 
calculated  immediately.  Since  the  inverse  problem 
has  only  a  few  unknowns,  reconstruction  of  target 
moments  is  overdetermined  and  least-squares  solu¬ 
tions  can  be  adopted.  For  each  target,  the  monopole 
moment  represents  the  integral  absorption,  the  dipole 
moments  result  in  correction  of  the  center  position, 
and  the  quadrupole  moments  lead  to  estimation  of 
the  target  volume,  in  whidi  tibe  integral  absorption 
can  be  redistributed.  The  target  volume  is  estimated 
by  simple  matching  of  the  quadrupole  moments  with 
those  of  an  ellipsoid,  which  has  the  freedom  to  rotate 
about  the  center  in  any  direction  to  any  angle. 

Figure  1  shows  a  schematic  of  the  combined  probe 
for  experiments,  which  is  a  circular  plate  made  of  rigid 
absorbing  material.  One  or  two  culraid  absorbers  were 
placed  2.5  cm  deep  in  a  homogeneous  mediiun.  For 
phantom  experiments,  raw  data  were  acquired  with 
our  combing  imaging  system.  The  ultrasound  array 
was  translated  in  the  x  direction  to  yield  target  vol¬ 
ume  estimates.  We  used  0.6%  Intralipid  as  a  homoge¬ 
neous  ba(^:ground,  which  yielded  a  reduced  scattering 
coefficient  of  6  cm~^  and  an  absorption  coefficient  of 
0.015  cm”^.  The  embedded  targets  were  as  scatter¬ 
ing  as  the  background  but  were  more  absorbing.  For 
simulations,  the  simulated  measurement  data  were  the 
sum  of  forward  model  prediction  and  random  noises 
ranging  from  0.1%  to  1%  (with  respect  to  the  ampli¬ 
tudes  of  the  incident  waves).  The  optical  properties 
of  homogeneous  media  and  targets  are  similar  to  those 
iised  in  phantom  experiments. 

Shown  in  Fig.  2  are  simulation  results  for  a  single 
target.  The  target  was  a  cuboid  of  (1.2,  0.8,  1)  cm 
in  the  x,  y,  and  z  directions,  respectively,  and  was 
centered  at  (0.2,  0,  2.5)  cm.  The  volume  of  the  target 
was  0.96  cm®.  The  absorption  coefficient  above  the 
background  was  0.25  cm“^.  Here  we  provide  the 
results  for  one  dipole  moment  (Dx)  and  one  quadruple 
moment  (Qxx)>  but  similar  results  were  obtained  for 
other  components.  The  initial  target  position  that 
was  input  into  the  reconstruction  algorithm  was  (0, 


0,  2.5)  cm,  which  deviated  slightly  from  the  true 
position.  This  initial  location  error  resulted  in  a 
negligible  effect  on  the  reconstructed  monopole  values, 
as  shown  in  Fig.  2(a).  The  offset  of  center  position 
can  be  corrected  by  the  estimated  dipole  moment, 
Dx  [see  Fig.  2(c)1,  which  is  slightly  lower  than  the 
true  value.  However,  the  inaccuracy  of  the  initial 
center  position  obviously  affects  the  reconstructed 
quadrupole  moment.  The  circles  in  Fig.  2(e)  repre¬ 
sent  values  greater  than  the  true  values.  After  we 
correct  the  center  position  according  to  the  estimated 
dipole  moments,  the  reconstructed  values  of  Qxx,  plot¬ 
ted  as  asterisks  in  the  same  figure,  are  much  closer  to 
the  true  values.  Although  the  mean  values  of  recon¬ 
structed  moments  are  essentially  independent  of  noise 


Fig.  1.  (a)  Bottom  view  of  the  combined  probe.  The 
central  rectangular  slot  gives  the  ultrasound  array  access 
to  tissues  underneath  the  probe.  The  circular  holes 
are  used  to  hold  optical  fibers.  The  small  holes  are  for 
light  sources,  and  the  larger  ones  are  for  detectors.  The 
diameter  of  the  probe  is  10  cm.  (b)  Side  view  of  the  probe. 


Noise  Level  xto*  Noise  Level  xto* 


Fig.  2.  Simulation  results  for  a  single  target,  (a),  (c),  and 
(e)  mean  values  of  monopole,  dipole,  and  quadrupole  mo¬ 
ments,  respectively,  (b),  (d),  and  (0  standard  deviations 
of  the  correspondi^  quantities.  The  solid  lines  in  (a),  (c), 
and  (e)  are  true  values. 
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Table  1,  Reconstructed  Absorption  Coefficients  from  Phnntnm  Experiments" 


Combination 

1 

Combination  2 

Value  (cm  ') 

Target  1 

Target  2 

Target  1 

Target  2 

True 

Reconstructed 

0.065 

0.064  ±  0.028 

0.065 

0.064  ±  0.027 

0.102 

0.106  +  0.026 

0.065 

0.064  ±  0.017 

"The  average  values  and  standard  deviations  are  estimated  finm  six  sets  of  experimental  data. 


level,  the  standard  deviations  of  all  moments  increase 
linearly  with  increasing  noise  level.  Higher-order 
moments  suffer  more  from  noise,  as  expected  intui¬ 
tively.  For  the  quadrupole  moments,  as  small  as  a 
0.2%  noise  level  could  lead  to  a  relative  error  of  ~20%. 
This  is  the  main  reason  that  we  did  not  attempt  to 
reconstruct  moments  beyond  the  second  order. 

To  test  the  capability  of  our  algorithm  to  charac¬ 
terize  multiple  targets,  we  performed  simulations  with 
two  absorbers  of  different  absorption  coefficients.  The 
absorbers  were  the  same  size  as  the  single  target.  One 
absorber  was  located  at  (—1.2,0, 2.5)  cm,  with  an  ab¬ 
sorption  coefficient  0.25  cm~^  beyond  the  backgroimd. 
The  other  absorber  was  located  at  (1.2, 0, 2.5)  cm,  with 
an  absorption  coefficient  0.1  cm~^  beyond  the  back¬ 
ground.  The  main  effect  of  adding  one  target  was  that 
the  relative  error  in  the  reconstructed  monopole  mo¬ 
ments  increased  to  ~2.5  times.  Less  of  an  effect  on  the 
dipole  moments  was  observed,  although  an  even  more 
trivial  difference  for  the  quadrapole  could  be  seen.  It 
was  also  found  that  the  noise-induced  deviation  was 
independent  of  target  optical  properties. 

Our  current  NIR  imaging  system  has  a  noise  level 
near  1.5%.  According  to  the  simulation  results,  the 
system  caimot  readily  retrieve  the  acciu-ate  quadrupole 
moments  from  measured  photon  density  waves.  Thus 
we  have  to  use  the  volume  information  obtained  from 
co-registered  ultrasound  images  in  target  characteri¬ 
zation.  We  used  two  geometrically  identical  targets 
that  were  approximately  1-cm*  cubes.  They  were  gel 
phantoms  made  from  0.6%  Intralipid  solution,  ink, 
and  ultrasound  scatterers.  The  boundaries  of  these 
targets  can  clearly  be  seen  in  ultrasound  images,  and 
target  centers  and  volumes  can  be  estimated  accu¬ 
rately.  For  example,  a  1.07-cm’  target  was  measured 
as  1.05  cm®  by  ultrasound.  For  co-registered  NIR 
imaging,  two  different  combinations  of  targets  were 
adopted.  The  first  was  a  pair  of  equally  absorbing 
cub^  with  an  absorption  coefficient  of  0.065  cm~^. 
The  second  was  a  pair  of  absorbing  cubes  with  absorp¬ 
tion  coefficients  of  0.065  and  0.1  cm“^,  respectively. 
We  simply  distribute  the  estimated  monopole  evenly 
over  the  target  volume.  So  the  reconstructed  ab¬ 
sorption  coefficient  was  an  average  value  within  each 
target.  The  true  and  the  reconstructed  vedues  are 
compared  in  Table  1.  The  mean  values  and  standard 
deviations  were  based  on  six  sets  of  measurements. 

Since  target  monopole  moments  can  be  estimated 
with  sufficient  accuracy,  the  reconstructed  absorption 
coefficient  depends  more  on  the  estimation  of  target 
volume,  which  can  be  obtained  from  the  reconstructed 
quadrupole  moment  or  co-registered  ultrasound  im¬ 
ages.  It  is  possible  that  the  extent  of  the  optical 
volume  of  a  lesion  is  larger  than  the  volume  esti¬ 


mated  from  the  acoustic  images.  This  is  one  of  om- 
motivations  for  estimating  the  optical  voliune  from 
the  quadrupole  moments.  Further  improvements 
are  needed  for  an  accurate  estimate  of  quadrupole 
moments.  In  the  moment-based  method  we  assume 
that  the  lesions  are  isolated  and  are  embedded  in  a 
homogeneous  background.  This  assumption  is  quite 
true  for  more-homogeneous  fatly  breasts  and  may  not 
hold  for  dense  breasts  that  consist  of  both  glandular 
tissue  and  fat.  In  the  latter  case,  we  could  segment 
ultrasound  images,  identify  tissue  typeSi  and  estimate 
background  optical  properties  of  different  tissues  in 
the  reconstruction  as  well.  Both  acoustic  and  optical 
contrasts  exist  in  tumors  but  the  sensitivities  of  these 
two  modalities  may  be  different.  Therefore,  the  cor¬ 
relation  between  acoustic  and  optical  heterogeneities 
remains  unknown  at  the  ciurent  stage  and  will  be 
determined  in  our  friture  clinical  studies. 

In  comparison  with  current  sophisticated  NIR  image 
reconstruction  algorithms,  our  method  is  very  simple 
in  terms  of  computation.  Once  an  ultrasound  image 
has  been  taken  and  targets  have  been  identified,  it 
is  possible  td  monitor  the  local  change  of  absorption 
continuously  and  in  real  time.  Averaging  over  a  large 
number  of  samples  will  help  suppress  the  effect  of  ran¬ 
dom  noise. 

In  this  Letter  an  algorithm  for  characterizing 
absorbers  has  been  demonstrated.  Theoretically,  a 
similar  algorithm  is  possible  for  characterizing  scat¬ 
tering  targets  as  well.  However,  since  the  weights 
of  monopole,  dipole,  and  quadrupole  moments  of 
scattering  targets  are  related  to  spatial  derivatives, 
which  are  two  order  higher  than  those  of  absorbers, 
the  quadrupole  moments  will  be  even  harder  to 
reconstruct. 
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Simultaneous  near-infrared  diffusive  light  and 
ultrasound  imaging 


Nan  Guang  Chen,  Puyun  Guo,  Shikui  Yan,  Oaqing  Piao,  and  Quing  Zhu 


We  have  constructed  a  near-real-tiiae  combined  imager  suitable  ibr  simultaneous  ultrasound  and  near- 
infirared  diffusive  li^  imaging  snd  coregistration.  The  imager  consists  of  a  combined  hand-held  probe 
and  the  associated  electroaka  far  data  acquisition.  A  two-dimensional  ultrasound  array  is  deployed  at 
the  center  of  the  combing  probe,  and  12  dual-wavelength  laser  source  fibers  (780  and  830  nm)  and  8 
optical  detector  fibers  are  deployed  at  the  periphery.  We  have  esperimentsMy  evaluated  the  efibeta  of 
missing  «qitical  sources  in  the  middle  rfthe  combined  probe  on  the  accuracy  ofthe  reconstructed  optical 
absorption  coeflBcient  and  anncissfid  the  improvements  of  a  reconstructed  absorption  coefficient  wi&  the 
guidance  of  the  coregistered  ultrasound.  ^Hie  results  have  shown  that,  wbrni  the  central  ultrasound 
array  area  is  in  the  nm^bboihood  2  cm  x  2  cm,  whidi  corresponds  to  the  size  of  most  conunercial 
ultrasound  transducers,  the  optical  imaging  is  not  affected.  The  results  have  also  shown  that  the 
iterative  inversion  algorithm  converges  quidcly  with  the  guidance  of  a  priori  three-dimensicmal  target 
distribution,  and  only  one  iteration  is  needed  to  reconstruct  an  accurate  optical  absorption  coefficient. 
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1.  Introduction 

Ultrasound  is  used  extensively  for  differentiation  of 
cysts  ffom  solid  lesions  in  breast  examinations,  and  it 
is  routinely  used  in  coiyunction  with  mammography. 
Ultrasound  can  detect  breast  l^ons  a  few  millime¬ 
ters  in  size.^  However,  its  specificity  in  breast  can¬ 
cer  diagnosis  is  not  considered  to  be  Ifi^  enough  as  a 
result  of  overlapping  characteristics  of  benign  and 
maliffoant  lesions.*^  Optical  imaging  based  on  dif¬ 
fusive  near-infiiared  (NIR)  lifi^t  has  the  great  poten¬ 
tial  to  diCEerentiate  tumors  finm  normal  breast 
tissues  throu^  determination  of  tissue  parameters, 
such  as  blood  volume,  blood  Og  saturatioit,  tissue 
light  scattering,  water  concentration,  and  the  concen¬ 
tration  and  lifetime  of  exogenous  contrast  agents.^-^^ 
As  a  potential  diagnostic  tool,  however,  NIR  difiusive 
light  imaging  suffers  firom  low  spatial  resolution  and 
lesion  location  uncertainties  because  of  intense  light 
scattering  in  tissue. 
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Most  NIR  imaging  reconstruction  algorithms  are 
based  on  tomographic  inversion  tediniques.^“ 
Reconstruction  of  tissue  optical  properties  in  general 
is  underdetermined  and  iU-pos^  because  the  total 
number  of  unknown  optical  properties  always  ex¬ 
ceeds  the  number  of  measurements,  and  the  pertur¬ 
bations  produced  by  the  heterogeneities  are  much 
smaller  than  the  badcground  sigoals.  In  addition, 
the  inversion  reconstruction  algorithms  are  sensitive 
to  measurement  noise  and  model  errors. 

Our  group  and  others  have  introduced  a  novel  hy¬ 
brid  imaging  method  that  combines  the  complemen¬ 
tary  features  of  ultrasound  and  NIR  difiusive  light 
imaging.*^"**  The  hybrid  imaging  obtains  coregis¬ 
tered  ultrasound  and  NIR  difbisive  li^t  images 
through  simultaneous  deployment  of  an  ultrasoimd 
array  and  NIR  source-detector  fibers  on  the  same 
probe.**-**-*^  Coregistration  permits  joint  evalua¬ 
tion  of  acoustic  and  optical  properties  of  breast  le¬ 
sions  and  enables  use  of  lesion  morphology  provided 
by  hig^-resolution  ultrasound  to  improve  the  lesion 
optical  property  estimate.  With  the  a  priori  knowl¬ 
edge  of  lesion  location  and  shape  provided  by  coreg¬ 
istered  ultrasound,  NIR  imaging  reconstruction  can 
be  localized  within  specified  three-dimensional  (3-D) 
regions.  As  a  result,  the  reconstruction  is  overdeter- 
minfld  because  the  total  number  of  unknown  optical 
properties  is  reduced  significantly.  In  addition,  the 
reconstruction  is  less  sensitive  to  noise  because  the 


1  December  2001  /  Vot.  40.  No.  34  /  APPUED  OPTICS  6367 


convergence  can  be  achieved  with  a  amnll  number  of 
iterations. 

The  clinical  use  of  the  combined  difl£pnn.ais  relies  on 
the  coregistration  of  both  ultrasound  and  NIR  sen¬ 
sors  at  the  probe  level.  Conventional  ultrasoimd 
pulse-echo  imaging  requires  that  an  imaging  trans¬ 
ducer  be  locat^  on  top  of  the  target,  whereas  NIR 
difiusive  light  imaging  is  feasible  when  the  optical 
source  and  detector  fibers  are  distributed  at  the  pe¬ 
riphery  of  the  ultrasound  tremsducer.  However,  the 
effects  of  missing  optical  sources  in  the  middle  of  the 
combined  probe  on  the  accuracy  of  the  reconstructed 
optical  properties  have  to  be  evaluated.  In  addition, 
the  improvements  of  reconstructed  optical  properties 
with  the  guidance  of  the  coregistered  udtrasound 
need  to  be  quantitatively  assessed.  Furthermore, 
real-time  data  acquisition  is  necessary  to  avoid  errors 
in  coregistration  caused  by  patient  motion  during  the 
dinifal  experiments.  In  this  paper  we  report  our 
experimental  results  on  the  optimal  probe  configura¬ 
tion,  and  we  quantify  the  improvements  on  recon¬ 
structed  optical  properties  using  a  combined  probe. 
We  also  demonstrate  simvdtaneous  combined  imag¬ 
ing  with  a  near-real-time  imager. 

2.  Near-Infrared  Diffusive  Wave  Imaging 

We  used  the  Bom  approximation  to  relate  the  scat¬ 
tered  field  I/ae'  (r,  <a)  measured  at  the  probe  surface  to 
absorption  variations  in  each  volume  element  within 
the  sample.  In  the  Bom  approximation,  the  scat¬ 
tered  wave  originated  fimm  a  source  at  r,i,  and  mea- 
stired  at  it  can  be  related  to  the  medium 
absorption  heterogeneity  Ajx.a(r„j)  at  Tyy  by 

[^tdluxi  ~  (1) 

where  Af  is  the  total  number  of  source-detector  pairs, 
N  is  the  total  number  of  imaging  voxels,  and  W^-  = 
G(r^;  r^i,  fci)I/iac(rv7.  w)vAryV^?  is  the  weig^it  ma¬ 
trix  given  in  Ref.  19.  G(r^j,  r^i,  w)  and  r^,  cj) 

are  a  Green’s  function  and  incident  wave,  respec¬ 
tively.  &>  is  the  modulation  frequency  and  D  is  the 
average  or  background  diffusion  coefficient,  which  is 
fhe  average  value  over  the  background  or  whole  tis¬ 
sue. 

With  M  measurements  obtained  fiiom  all  possible 
source-detector  pairs  in  the  planar  array,  we  can 
solve  N  unknowns  of  by  inverting  the  above  matrix 
equation.  In  general,  the  perturbation  Eq.  (1)  is  un¬ 
derdetermined  (M  <  N)  and  ill-posed. 

NIR  imsiging  by  itself  generally  has  poor  depth 
discrimination.  However,  iiltrasound  is  excellent  in 
providing  accurate  target  depth.  Once  the  target 
depth  is  available  firom  coregistered  ultrasound,  we 
can  set  of  a  nontarget  depth  equal  to  zero.  This 

implies  that  all  the  measured  pertt^ations  originate 
firom  the  particular  depth  that  contains  the  target. 
Because  the  number  of  unknowns  is  reduced  signifi¬ 
cantly,  the  reconstaruction  converges  very  fast.  In 
Ref.  23  we  reported  that,  with  a  priori  target  depth 
provided  by  ultrasound,  the  accuracy  of  the  recon- 
stmcted  has  been  improved  by  15-30%  on  aver- 
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age,  and  the  speed  of  reconstmction  has  been 
improved  by  an  order  of  magnitude.  In  this  paper 
we  furthermore  demonstrate  that,  with  the  3-D  tar¬ 
get  distribution  provided  by  coregistered  ultrasound, 
the  accuracy  of  reconstmcted  p.,,  and  the  reconstmc¬ 
tion  speed  can  be  further  improved. 

To  solve  the  imknown  optical  properties  of  Eq.  (1), 
we  used  the  total  least-squares  (TUS)  method^-*^  to 
iteratively  invert  Eq.  (1).  The  TLS  method  performs 
better  than  other  least-squares  when  the  measime- 
ment  data  are  subject  to  noise  and  the  linear  operator 
W  contains  errors.  We  found  that  the  TLS  method 
provides  more  accurate  reconstmcted  optical  proper¬ 
ties  than  other  least-squares  methods,  and  we 
adopted  TLS  in  solving  inverse  problems.  It  has 
been  shown  in  Ref.  28  that  the  TLS  minimization  is 
equivalent  to  the  following  minimization  problem: 


11x11* +1  ’ 


where  X  represents  unknown  optical  properties. 
The  conjugate  gradient  technique  was  employed  to 
iteratively  solve  Eq.  (2). 

3.  Methods 


A.  Combined  Probe  and  Imaging  Geometry 
There  are  four  basic  requirements  to  guide  the  design 
of  the  combined  probe.  First,  reflection  geometry  is 
preferred  because  a  conventional  ultrasound  scan  is 
performed  with  this  geometry.  Second,  an  ultra¬ 
sound  array  needs  to  occupy  tiie  center  of  the  com¬ 
bined  prol^  for  coherent  imaging.  Third,  NIR 
sources  and  detectors  have  to  be  distributed  at  the 
periphery.  Because  photon  propagation  distribu¬ 
tion  exhibits  a  banana  shape,  imaging  of  the  tissue 
volume  underneath  the  probe  is  feasible  even 
through  there  are  no  sources  and  detectors  deployed 
in  the  central  portion  of  the  probe.  Fourth,  the  min¬ 
imum  source-detector  separation  should  be  larger 
than  1  cm  for  the  diffusion  approximation  to  be  valid, 
and  the  maximum  separation  should  be  ~8-9  cm  to 
effectively  probe  deptiis  of  3-4  cm. 

On  the  basis  of  these  requirements  we  deployed  12 
dual-wavelength  optical  source  fibers  and  8  detector 
fibers  over  a  9  cm  X  9  cm  probe  area  (see  Fig.  1). 
The  tniniTniiTn  and  maximum  source-detector  sepa¬ 
rations  in  the  configuration  are  1.4  and  8  cm,  respec¬ 
tively.  To  study  the  effect  of  the  centred  optical  hole 
on  the  accuracy  of  the  reconstructed  optical  proper¬ 
ties,  we  compared  the  reconstruction  results  with  an 
extra  center  source  and  without  the  center  source. 
The  configuration  without  the  center  source  corre¬ 
sponds  to  a  2  cm  X  2  cm  hole  area.  We  further 
moved  the  noncenter  12  sources  and  8  detectors  to¬ 
ward  periphery  by  leaving  a  3  cm  x  3  cm  hole  area  in 
the  middle.  Figure  2  shows  the  picture  of  a  com¬ 
bined  probe  with  the  3  cm  X  3  cm  central  area  occu¬ 
pied  by  an  ultrasoimd  array.  The  ultrasound  array 
consists  of  64  elements  made  of  1.5-mm-diameter  pi¬ 
ezoelectric  transducers  (Valpey  Fisher  Inc).  The 


*  H  "  fcm  ~  '  >1 

Fig.  1.  Schematic  airangement  of  NIR  source  and  detector  fibers 
on  the  |Mt)be.  Small  solid  circles  are  the  source  fibers  and  larger 
■olid  ejdes  are  the  detector  fibers. 


transducers  are  deployed  in  a  rectangular  matrix 
xrith  4-nim  spacing  in  both  x  and  y  dir^:tions.  The 
center  firequency  of  the  transducer  is  6  MHz  and  the 
bandwidth  is  40%.  The  transducers  are  made  firom 
the  same  piece  of  piezoelectric  transducer  material, 
llierefore  the  gain  difference  among  different  trans¬ 
ducers  is  less  than  3  dB.  The  12  dual-wavelength 
optical  laser  diode  sources  (760  and  830  nm)  wd  8 
photomultiplier  tube  (PMD  detectors  are  coupled  to 
.the  probe  through  optical  fibers,  which  are  deployed 
at  the  periphery  of  the  two-dimensional  (2-D)  ultra¬ 
sound  array.  This  hybrid  array  deployment  compro¬ 
mises  ultrasound  coherent  imaging  and  NIR  diffusive 
li^t  imaging  characteristics. 

The  9cmX9cmX4cm  image  volume  underneath 
the  probe  is  discretized  into  voxels  of  size  0.4  cm  X  0.4 
cm  X  1  cm.  There  is  a  trade-off  between  the  accu¬ 
rate  estimation  of  the  weight  matrix  W  and  the  voxel 


Fig.  2.  Picture  of  an  ea^rim'ental  probe.  Aniiltrasouiuiarrayof 
8  X  8  »  64  transducers  occupies  the  central  3  cm  X  3  cm  area,  and 
12  dual-wavelength  source  fibers  and  8  detector  fibers  are  deployed 
at  the  periphery.  •' 


Fig.  3.  Schematic  of  the  NIR  firequenqr-domain  imaging  system. 
The  modulation  firequengr  is  140  MHz.  The  12  dual-wavdength 
source  channels  are  switched  on  sequentially  by  a  PC,  and  8  de¬ 
tector  channels  receive  signals  in  paraUeL  BPF,  bandpass  filter; 
OSC,  oscillator. 

size.  Because  Wy  is  a  discrete  approximation  of  the 
integral 

J  G(r„  Trf,  o))t/fac(r„  r„  w)  ^  dr?, 

it  is  more  accurate  when  ffie  voxel  size  is  smaller. 
However,  the  total  munber  of  reconstructed  un¬ 
knowns  will  increase  dramatically  wiffi  the  decreas¬ 
ing  voxel  size.  Furthermore,  the  rank  of  W  does  not 
increase  in  the  same  order  as  .the  total  number  of 
voxels  when  the  voxel  size  decreases.  This  suggests 
that  neighboring  W^s  are  correlated  when  the  voxel 
size  is  smaller,  and  a  further  decrease  in  voxel  size 
will  not  add  more  independent  information  to  the 
weight  matrix.  We  found  that  a  0.4  cm  X  0.4  cm  x 
1  cm  voxel  size  is  a  good  compromise.  Therefore  we 
used  this  voxel  size  in  image  reconstructions  reported 
in  this  paper. 

B.  Experimental  Systems 

1.  Near-Infrared  Imaging  System 
We  constructed  a  NIR  fi%quency-domaia  imaging 
system.  The  block  diagram  of  the  system  is  shown 
in  Fig.  3.  This  system  has  12  dual-wavelength 
source  chaimels  and  8  pauallel  receiving  channels. 
On  the  transmission  part,  12  pairs  of  dual¬ 
wavelength  (780  and  830  nm)  laser  diodes  are  used  as 
light  sources,  and  their  outputs  are  amplitude  mod- 
u^ted  at  140.000  MHz.  Eadi  one  of  the  12  optical 
combiners  (OZ  Optics  Inc.)  looks  like  a  Y  adapter, 
guiding  the  emission  of  two  diodes  of  different  wave- 
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lengths  through  the  same  thin  optical  fiber  (approx¬ 
imately  0.2  mm  in  diameter).  To  reduce  noise  and 
interference,  an  individual  driving  circuit  is  built  for 
each  diode.  As  a  laser  diode  works  in  series,  a  con¬ 
trol  board  that  inteiprets  instructions  fi:um  a  PC  is 
used  to  coordinate  operations  of  associated  compo¬ 
nents.  When  a  sin^e  transmission  chaTini*!  is  se¬ 
lected,  it  turns  on  the  corresponding  driving  circuit  so 
that  a  dc  driving  current  can  be  set  up  for  the  diode. 
At  the  same  time,  a  selected  signal  is  sent  to  a  rf 
switching  unit,  wWch  distributes  a  rf  signal  to  the 
ri|^t  channel  to  modulate  the  optical  output.  On 
the  reception  part,  eight  PMTs  are  employed  to  detect 
diffusely  reflected  light  fi:om  turbid  media.  Each 
PMT  is  housed  in  a  sealed  aluminum  box,  shielding 
both  environment  lights  and  electromagnetic  fields, 
and  an  optical  fiber  (3  mm  in  diameter)  couples  NIR 
light  from  the  detection  point  to  the  reception  window 
of  the  PMT.  The  electrical  signal  converted  firom  the 
optical  input  is  generally  weak  and  rather  hi^  in 
fi^uen<7,  so  hi^-gain  amplification  and  firequency 
transform  are  necessary  before  it  can  be  sampled  by 
an  anal(%-to-digital  (A^)  board  inside  the  PC.  We 
built  ei^t  parallel  heterodyne  amplification  dian- 
nels  to  measiue  the  response  of  all  detectors  simul¬ 
taneously,  which  reduces  the  data-acquisition  time. 
Each  amplification  channel  consists  of  a  rf  amplifier 
(40  dB),  a  mixer  in  which  the  rf  signal  (OSCl,  140.000 
MHz)  is  mixed  with  a  local  oscillator  (OSC2, 140.020 
MHz),  a  bandpass  filter  centered  at  20  kHz,  and  a 
low-fi:equen<7  amplifier  of  30  dB.  The  heterodyned 
two-stage  amplification  scheme  helps  suppress  wide¬ 
band  noises  efficiently.  We  also  generated  a  refer¬ 
ence  signal  of  20  kHz  by  directly  mixing  OSCl  and 
OSC2,  which  is  necessary  for  retrieving  phase  shifts. 
Ei|^t  detection  signals  and  one  reference  are  sam¬ 
ple  converted,  and  acquired  into  the  PC  simulta¬ 
neously,  in  which  the  Hilbert  transform  is  used  to 
compute  the  amplitude  and  phase  of  each  chaimel. 
The  entire  data  acquisition  takes  less  than  1  min, 
which,  is  fast  enough  to  acquire  data  fipom  patients. 

One  of  the  challenges  encountered  in  the  design  of 
a  NIR  imaging  system  is  the  huge  dynamic  range  of 
signals  received  at  various  source-detector  dis¬ 
tances.  For  example,  for  a  semi-infinite  phantom 
made  of  0.5%  Intralipid  solution,  the  amplitude  mea¬ 
sured  at  1  cm  away  firom  a  source  is  approximately 
5000  times  larger  ffian  that  at  8-cm  separation. .  In 
addition,  the  perturbation  that  is  due  to  an  embedded 
heterogeneity  with  optical  properties  similar  to  a  tu¬ 
mor  is  normally  a  few  percent  of  the  background 
Bignfll-  As  a  result,  a  reflection-mode  NIR  imaiging 
system  should  have  at  least  a  120-dB  dynamic  range 
to  probe  a  target  up  to  4  cm  in  depth.  It  is  hard  to 
build  amplifiers  that  work  linearly  over  such  a  wide 
dynamic  range.  We  overcome  this  difficulty  by  im¬ 
plementing  two-level  source  outputs.  The  dc  output 
of  a  laser  diode  is  controlled  when  its  feedback  loop  is 
a<liusted,  whereas  the  rf  signal  is  switched  simulta¬ 
neously  by  a  two-step  attenuator  (no  attenuation  or 
30-dB  attenuation).  When  the  source  and  detector 
are  dose  to  each  other,  the  source  is  controlled  to 


yield  a  low-level  output.  When  the  separation  be¬ 
comes  larger,  a  30-dB  higher  output  level  shoidd  be 
used.  With  this  two-level  source  scheme,  our  system 
achieved  fairly  good  linearity  over  a  wide  range  of 
source-detector  separations  (fix)m  1.5  to  8  cm). 

Because  the  parameters  of  an  individual  laser  di¬ 
ode  or  a  PMT  vary  considerably  firom  one  to  another, 
we  have  to  calibrate  tiie  gain  and  phase  shift  for  each 
channel.  A  set  of  measurements  obtained  firom  all 
source-detector  pairs  placed  on  the  boundary  of  a 
homogeneous  medium  is 

Aap,  a  =  1,  2 . m;  p  =  1,  2,  . . . ,  n. 

Here,  ampUtude  and  phase  are  related  to 
source  a  and  detector  p,  and  m  and  n  are  the  total 
number  of  sources  and  detectors,  respectively. 
From  the  diffusion  theory,  we  can  obtain  the  follow¬ 
ing  set  of  equations^: 

A  —Tt  NT  /r»\ 

Pop 

“  <Pi(«)  +  ViiO)  +  *rP«p, 

in  which  I,{a)  and  <p,(a)  are  the  relative  gain  and 
phase  delay  associate  with  soiurce  chaimel  a,  IJff) 
and  (p,f(P)  are  similar  quantities  associated  with  de¬ 
tector  channel  p,  is  the  corresponding  separation, 
and  kr  +  jki  is  tiie  complex  wave  number.  We  obtain 
the  following  set  of  linear  equations  by  taking  a  log¬ 
arithm  of  the  above  equations  related  to  amplitude: 

log(p„p*A„p)  «  log{/,(a)]  +  log[/rf(p)]  -  kiP^, 

$ap  =  q>,(a)  +  <PclO)  +  *rP«p-  (3) 

Although  the  optical  properties  of  the  calibration  me¬ 
dium  are  known  in  advance,  we  leave  the  wave  num¬ 
ber  as  a  variable  and  use  fitted  kr  and  ki  to  calculate 
the  background  scattering  and  absorption  coeffi¬ 
cients.  We  verified  our  calibration  meffiod  by  com¬ 
paring  the  best  fitted  k’s  with  real  values.  The 
results  of  our  using  0.5-0.8%  Intralipid  solutions  al¬ 
ways  yielded  scattering  and  absorption  coefficients 
with  a  rather  good  accuracy.  With  the  two  unknown 
wave  numbem  included,  the  total  number  of  un¬ 
knowns  is  2(in  ->*11  +  2),  which  is  generally  far 
smaller  than  the  number  of  measurements  m  X  n. 
Consequently,  Elq.  (3)  is  overdetermined.  We  can 
solve  all  /,(“)» <PdO)  terms  as  well  as 
two  unknown  wave  numbers  in  a  least-squares  sense. 
Then  all  measurements  can  be  calibrated  accord¬ 
ingly.  The  results  of  amplitude  A^p  =  exp(-A!iP„p)/ 
Pap^  and  phase  <|>„p  =  i^Pap  after  calibration  are 
shown  in  Fig.  4.  As  one  can  see,  the  calibrated  am¬ 
plitude  Gog  Pap*A„p)  and  the  phase  firom  various 
source-detector  pairs  change  linearly  with  distance. 

2.  Ultrasound  System 

The»ultrasound  system  diagram  is  shown  in  Fig.  5, 
and  the  system,  consists  of  64  parallel  transmission 
and  receiving  chaimels.  Each  transmission  circuit 
dan  generate  a  high-voltage  pulse  of  200-ns  duration 
(6  bniz)  with  125  V  peak  to  peak  to  the  connected 
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(b) 

FSg.  4.  (a)  Log(Paa*^«a)  versus  distancev  after  calibratioa. 
(b)  Phase  versus  distance  after  calibration. 


transducer.  Each  receiving  circuit  two-stage 
amplifiers  followed  by  an  A/D  converter  with  40-M]^ 
sampling  firequency.  -The  amplifier  gain  can  be  con¬ 
trolled  based  on  the  target  streng^.  A  group  of 
transmission  channels  can  be  addressed  simulta¬ 
neously  to  transmit  pulses  firom  neighbor  transduc¬ 
ers  with  specified  delays  and  therefore  to  focus  the 
transmission  beam.  The  retuned  signals  can  be  re¬ 
ceived  simultaneously  by  a  group  of  transducers,  and 
the  signfli.s  can  be  sununed  with  specified  delays  to 
form  a  receiving  beam.  » ■  •  .  « •  ; 

The  data-acquisition  procedure  is  the  following. 
The  first  9-element  neighbor  subarray  (dashed- rect¬ 
angle  in  jFig.  6)  firom  the  64-element  transducer  array 
flnH  the  corresponding  channels  are  chosen^ andthen 
the  transmission  delay  profiles  are  generated  in  the 
computer  according  to  the  prespecified  focal  depth. 
The  delay  profile  data  are  transferred  to.  the:, 64- 
channel  delay  profile  generator,  which  tri^ers  the64 
high-voltage  pulsers  as  well  as-  the  re^ivingi  chan¬ 
nels.  The  returned  ultrasound  signals  are  amplified 


-Tl 


Fig.  5.  Schematic  of  our  ultrasound  scanner.  We  connected  64 
ultrasound  transducers  to  64  parallel  transmission  and  reception 
channels.  The  transmission  part  consists  of  64  high-voltage 
pulsers,  which  can  be  controlled  by  computer-generated  delay  pro¬ 
files.  The  reception  part  consists  of  64  two-stage  amplifiers  and 
A/D  converters.  CH,  channel 


by  two-stage  amplifiers  and  sampled  by  A/D  convert¬ 
ers.  The  data  are  buffered  in  the  memories  and  are 
read  by  the  computer  after  the  entire  data  acquisition 
is  completed.  The  second  subarray  (solid  rectangle 
in  Fig.  6)  is  chosen  and  the  same  data-acquisition 
process  is  repeated.  A  total  of  64  subarrays  is  used 
in  the  data  acquisition.  After  the  64-subarray  data 
acquisition  is  completed,  the  data  stored  in  the  mem¬ 
ories  are  read  by  the  computer  for  image  formation. 
The  entire  data  acquisition  and  imaging  display  are 
performed  in  approximately  5  s,  whi<^  is  fast  enough 
for  clinical  erqperiments.  To  ensure  good  signal-to- 
noise  ratio,  we  perform  all  the  electronics  using 
printed  circuit  boards. 

Figure  7  shows  the  pictuio  of  the  entire  system  and 
the  combined  probe.  Both  the  NIR  system  (top)  and 
the  ultrasound  system  (bottom)  are  mounted  on  a 
hospital  cart.  The  combined  probe,  which  houses 
the  ultrasound  array  and  the  NIR  source-detector 
fibers,  is  designed  to  be  hand  held  to  scan  patients. 
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Fig.  6.  Ultrasound  subarray  scanning  configuration. 
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Fig.  7.  Picture  of  our  combined  system.  NIR  system  (top)  and 
ultrasound  ^tem  (bottom)  are  mounted  on  a  hospital  cart. 


C.  Phantoms 

We  used  0.5-0.6%  Intralipid  solutions  to  mimic  nor¬ 
mal  human  breast  tissues  in  all  experiments,  and  the 
corresponding  reduced  scattering  coefiBcient  p.,' 
ranges  from  5  to  6  cm~^  The  Intralipid  is  contained 
in  a  large  fish  tank  to  set  up  approximately  a  semi¬ 
infinite  homogeneous  phantom.  Small  semispheri- 
cal  balls  (1  cm  in  diameter),  made  of  acrylamide  gel,“ 
are  inserted  into  Intralipid  to  emulate  lesions  emb^- 
ded  in  a  breast.  The  r^uced  scattering  coefficients 
of  the  gel  phantoms  are  similar  to  that  of  the  back¬ 
ground  m^um  (p,,'  *•  6  cm"^),  and  we  changed  the 
absorption  coefficients  to  difierent  values  by  adding 
different  concentrations  of  India  ink  to  emulate  high- 
contrast  (pa  =  0.25-cm“^)  and  low-contrast  (p^  = 
0.1-cm“*)  lesions.  Ultrasound  scattering  particles 
of  200  pm  in  diameter  are  added  to  the  gel  phantom 
before  the  gel  is  formed. 

4.  Experimental  Results 

A.  Effects  of  Missing  Optical  Sources  in  the  Combined 
Probe 

A  series  of  experiments  was  conducted  to  estimate 
the  optimal  hole  size.  Three  probe  configurations 
were  investigated:  (a)  no-hole,  (b)  2  cm  X  2  cm  cen¬ 
tral  hole,  and  (c)  3  cm  X  3  cm  hole  probes.  The 
no-hole  probe  was  essentially  the  same  as  case  (b) 
except  that  an  additional  light  source  was  added  in 
the  middle.  Figmre  8  shows  reconstructed  NIR  im¬ 
ages  for  on-center  targets  of  high  (p^  =  0.25  cm“S  left 
column)  and  low  contrast  (p^  =  0.1  cm“^,  right  col- 
tuxm)  located  2.5  cm  deep  inside  the  Intralipid.  The 
fitted  background  p„  and  p,  are  0.015  and  5.36  cm” 
respectively.  Witii  the  tiuget  depth  provided  by  ul¬ 
trasound,  we  performed  reconstruction  in  the  target 
layer.  The  centers  of  the  voxels  in  this  layer  were  (x, 
y,  2.5  cm),  where  x  and  y  were  discrete  spatial  x-y 
coordinates,  and  the  thic^ess  of  the  layer  was  1  cm. 
For  the  hi^-contrast  target  case,  there  are  no  im¬ 
portant  differences  in  image  qtiality  associated  with 
different  probes  [Figs.  8(a)  and  8(c)]  except  that  with 
a  3  cm  X  3  cm  hole.  The  first  row  of  Table  1  provides 


measiired  maximum  p„  values  fiem  the  correspond¬ 
ing  images.  Because  of  the  low  spatial  resolution  of 
diffusive  imaging,  the  iMundaries  of  the  targets  are 
not  well  defined.  The  maximum  value  is  a  better 
estimation  of  reconstructed  target  p„.  From  no  hole 
to  2  cm  X  2  cm,  the  reconstructed  maximum  p^  de¬ 
creases  slowly.  But  for  3  cm  X  3  cm,  the  maximiiTn 
p„  drops  suddenly  to  0.104  cm~^,  which  is  less  than 
half  of  the  origin^  value.  Another  imaging  param¬ 
eter  we  measured  is  the  full  width  at  half-maximiun 
(FWHM)  of  the  corresponding  images.  Because  the 
image  lobes  were  elliptical  in  general,  we  measured 
the  widths  of  longer  and  shorter  axes  and  used  the 
geometric  mean  to  estimate  the  FWHM.  The  re¬ 
sults  are  shown  in  Table  1,  and  the  FWHM  almost 
increases  with  the  hole  size.  We  also  measured  the 
image  artifact  level,  which  was  defined  as  the  ratio  of 
the  peak  artifact  to  the  maximum  strength  of  the 
image  lobe  and  is  given  in  decibels.  The  results  are 
shown  in  Table  1.  No  artifacts  were  observed  in  tiie 
images  of  no-hole  and  2  cm  x  2  cm  hole  probes. 
However,  the  peak  artifact  level  at  the  -  14.3-<ffi  level 
was  measured  in  the  image  of  the  3  cm  x  3  cm  hole 
probe.  When  the  contrast  was  low,  the  recon¬ 
structed  maximum  absorption  coefficients  and  mea¬ 
sured  FWHMs  were  essentially  the  same  for  the 
no-hole  and  2  cm  X  2  cm  hole  probes.  However,  the 
reconstructed  maximum  value  dropped  to  60%  of  the 
true  value  for  the  3  cm  x  3  cm  probe.  The  artifact 
levels  meastired  in  the  images  of  three  probe  config¬ 
urations  were  similar  and  were  worse  than  the  high- 
contrast  case.  The  image  artifacts  are  related  to  the 
reconstmction  algorithm.  When  the  target  contrast 
is  weak  or  the  signal-to-noise  ratio  is  low,  the  inver¬ 
sion  algorithm  pi^uces  artifacts  around  ffie  edges  of 
the  images. 

For  shallow  targets  (here  we  set  the  target  depth  to 
be  1.5  cm)  the  NIR  system  has  a  relatively  poorer 
performance.  This  is  due  to  less  source-detector 
pairs  experiendng  the  existence  of  a  shallow  absorber. 
As  shown  in  Fig.  9,  image  artifacts  are  obviously  worse 
compared  with  Fig.  8.  However,  the  conclusion  about 
the  hole  size  of  the  probe  remains  true.  Table  2  lists 
all  the  measur^  imaging  parameters  obtained  firom 
three  probe  configurations.  Although  a  3  cm  X  3  cm 
hole  is  somewhat  too  big  to  obtain  go^  enou^  results, 
the  optimal  hole  size  is  in  the  nei^bhoihood  of  2  cm  X 
2  cm.  This  optimal  size  is  approximately  the  size  of 
commercial  ultrasound  transducers. 

In  the  above  studies,  we  used  the  iteration  number 
obtained  from  the  no-hole  configuration  for  the  rest  of 
the  configurations.  Ideally,  the  iteration  should 
stop  when  the  object  function  [see  Eq.  (2)]  or  the  error 
peiformance  sxuface  reaches  the  noise  floor.  How¬ 
ever,  system  noise,  particularly  coherent  noise,  was 
difficult  to  estimate  fa)m  experimental  data.  In  gen¬ 
eral,  we  fotmd  that  the  reconstructed  values  were 
closer  to  true  values  when  the  object  function  reached 
approximately  5-15%  of  the  initial  value  (total  en¬ 
ergy  in  the  measurements).  Therefore  we  used  this 
criterion  (~1()%  of  the  initial  value)  for  ffie  no-hole 
configvuration.  Because  the  signal-to-noise  ratio  of 
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Fig.  8.  Reconstructed  NIR  images  of  deepw  ibLiget^(2^  depth,  1  cm  in  diameter,  and  the  fitted  background  iv.  and  p.'  are  0.015 
and  5.36  cm~\  respectively).  The  left  ooluinn.oonmpo^iib  tomiages  of  a  high-contrast  target  (|i.„  =  0.25  cm~^)  obtained  firom  different 
probeconfigurations,andtherig^tcolumncon;M]^^'to.ijma^ofaiow-contrasttarget(M„  •0.1cm~^).  Each  row  is  related  to  a  specific 
hole  size:  (a)  and  (b)  no  hole,  (c)  and  (d)  2  cm  X  (^)  and  (P  3  cm  X  3  cm. 


.  ■■ sfii  »\)Uy 

the  data  decreased  with  the  increase  in  hole  si^,  yre 
could  not  find  consistent  criterion  for  ^th  jid-hple 
and  hole  data.  Therefore  we  used  the  stime  iteration 
number  obtained  firom  the  no-hole  case  for  the^^pjb 
configurations,  and  the  comparison  was,  bpspdi»n|iie 
same  iteration  number.  ... 

•  ••  ,  .ras^Oti^q 

B.  Ultrasound-Guided  Near-Infrared  lnnagjng.'i|  .^f^n^ 

Three-dimensional  vdtrasoimd  images  can*.;pitOTi& 
3-D  distributions  of  targets.  With  the  dprioraatarget 


depth  information,  the  optical  reconstruction  can  be 
improved  significantly.  An  example  is  given  in  Fig. 
10.  The  target  again  was  a  1-cm-diameter  gel  ball  of 
low  (jia  .=  O.lrcm"^)  optical  contrast  and  was  embed¬ 
ded  :at  approximately  (0,  0,  2.5  cm)  inside  the  In¬ 
tralipid  medium.  The  fitted  background  p.a  and  p.,' 
•are  0.02  and  5.08  cm~\  respectively.  The  combing 
^robe  shown  in  Fig.  2  was  used  to  obtain  the  ultra¬ 
sound  and  NIR  data  simultaneously.  Figure  10(a) 
■ahows  an  A-scan  line  of  a  returned  ultrasound  echo 
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TaUal.  ParametefS  of  Raconstmctedlmagas  for  Oaap  High-Contrast 
(it,  -  0.25-cni-^)  and  LowContraat  (|«,  m  ai-cm"’)  Targata* 

_ Prolw'firpe 


Parameter 

No  Hole 

2  cm  X  2  cm 

3  cm  X  3  cm 

Hidi  contrast 

«  O 

0.251 

0.234 

0.104 

(cm"*) 

FWHM' 

(cm) 

1.85 

1.91 

2.44 

Artifacts 

(dB) 

Background 

(-22) 

Background 

-14.3 

Low  contrast 

iVuuc) 

(cm"*) 

0.105 

0.111 

0.064 

FWHM 

(cm) 

2.02 

1.83 

2.16 

Artifacts 

(dB) 

-6.90 

-8.10 

-5.65 

*n>e  fitted  background  |iu  and  |t,'  are  0.015  and  5.36  cin~S 
respectively. 

measured  maximum  ..value  of  the  reconstructed 
absorption  coefiSdent  map. 

*FWHM  is  defined  as  the  geometric  mean  of  the  widths  mea¬ 
sured  at  longer  and  shorter  axes  of  the  elliptical  image  lobe. 


signal  received  by  one  ultrasound  transducer  located 
on  top  of  the  target.  As  acoustic  scatters  were  uni¬ 
formly  distributed  in  the  target,  signals  were  re¬ 
flected  from  inside  the  target  as  well  as  from  the 
surfaces.  The  reflected  signals  from  the  front  and 
back  surfaces  of  the  gel  ball  can  be  clearly  identified 
in  the  echo'  signal.  On  the  basis  of  the  target  depth, 
we  reconstructed  the  optical  absorption  coefficient  at 
the  target  depth  only  (1  cm  in  thickness)  by  setting 
the  purtrubations  frpom  the  other  deptl^  equal  to 
zero.  We  also  performed  3-D  optic^-oifly  recon¬ 
struction.  Figure  10(b)  shows  the  reconstructed  ab¬ 
sorption  image  frrom  a  S-D  optical-only  reconstruction 
[layer  three  of  voxel  cpordi^tes  (x,  y,  2.5  cm)  and  1 
cm  thick],  whereas  Fig;  10(c)  8how;8  the  reconstructed 
image  of  the  corresponding  target  fruin  ultrasound- 
guided  reconstruction.  For  opti^-k>nly  reconstruc¬ 
tion,  the  algorithm  did  not  converge  to  a  localized 
spatial  r^on,  and  the  image  contrast  was  poor. 
*^0  measured  maximum  absorption  coefficient  was 
0.088  cm~\whidi  was  close  to  the  true  value.  How¬ 
ever,  the  measured  spatial  location  of  the  maximum 
value  was  (—1.6,  -1.2  cm),  which  was  too  far  frnm  the 
true  target  location.  With  the  a  priori  target  depth, 
the  reconstruction  performed  at  ffie  target  layer  can 
localize  the  target  to  the  correct  spatial  position. 
The  measured  maximum  absorption  coefficient  was 
0.12  cm“^  and  its  location  was  (0, 0.4  cm),  which  was 
very  dose  to  the  true  target  location.  example 
demonstrates  that  a  priori  target  depth  can  signifi¬ 
cantly  improve  the  reconstruction  accuraf^  and  tar¬ 
get  localization. 

In  addition  to  use  of  a  priori  target  depth  informa¬ 
tion,  we  can  also  use  the  target  spatial  distribution 
provided  by  ultrasound,  to  guide  ffie  reconstruction. 
We  performed  a  set  of  experiments  with  two  targets 


located  at  2.5  cm  in  depth  inside  the  Intralipid. 
Each  target  is  a  1-cm^  gel  cube  containing  ultrasound 
scatters.  For  optical  properties,  they  both  could  be 
high  contrast  ((!„  =  0.25  cm“^)  or  low  contrast  (p,„  = 
0.1  cm~^),  but  had  the  same  r^uced  scattering  coef- 
fident  as  the  background.  The  fitted  background  (Iq 
and  are  0.017  and  4.90  cm“\  respectively.  One 
target  was  centered  approximately  at  (- 1.0,  - 1.0, 2.5 
cm),  whereas  the  other  was  at  (1.0, 1.0, 2.5  cm).  The 
distance  between  the  centers  of  the  two  targets  was 
2.8  cm. 

Figure  11(a)  is  the  ultrasoimd  image  of  two  high- 
contrast  targets.  As  the  field  of  view  of  the  ultra¬ 
sound  system  was  nearly  a  3  cm  X  3  cm  square,  these 
two  targets  appeared  at  diagonal  comers.  The  mea¬ 
sured  peak  positions  of  the  two  targets  were  (-0.6, 
—1.0  cm)  and  (1.0,  1.0  cm),  which  differed  frrom  the 
true  target  locations  by  only  one  voxel.  The  low  con¬ 
trast  of  ffie  ultrasound  image  is  related  to  the  spedde 
noise.  Because  our  ultrasound  array  is  sparse,  the 
imaging  quality  is  not  state  of  the  art  (see  more  dis¬ 
cussion  in  Section  5).  The  NIR  image  of  these  tar¬ 
gets  was  obtained  simultaneously  and  is  shown  in 
Figure  11(b).  We  performed  the  reconstruction  at 
the  target  layer  by  taking  advantage  of  target  depth 
information.  A  total  of  123  iterations  was  used  to 
obtain  Fig.  11(b).  The  measured  peak  positions  of 
the  two  targets  were  (—1.4,  -1.0  cm)  and  (0.6,  0.6 
cm),  which  were  one  voxel  off  frrom  tiie  true  target 
‘  locations  (-1.0,  -1.0  cm)  and  (1.0,  1.0  cm),  respec¬ 
tively.  The  corresponding  reconstmcted  absorption 
coefficients  were  0.242  and  0.251  cm"^,  which  were 
dose  to  the  true  values.  However,  the  two  targets 
were  almost  connected  to  each  other,  and  their  spa¬ 
tial  localization  was  poor.  For  low-contrast  targets, 
the  ultrasoimd  image  is  shown  in  Fig.  11(c),  and  the 
measured  peak  locations  of  the  two  taints  were 
{— 1.0,  —0,6  cm)  and  (0.6, 1.0  cm),  which  differed  frrom 
the  true  target  locations  by  only  one  voxel.  The  cor¬ 
responding  NIR  image  is  shown  in  Fig.  11(d),  and  the 
measured  peak  locations  of  the  two  targets  were 
(-2.2,  —1.0  cm)  and  (0.6,  1.0  cm).  The  left  target 
was  off  the  true  location  by  three  voxels.  The  corre¬ 
sponding  reconstructed  absorption  coeffidents  were 
0.063  and  0.1004  cm~^  at  87  iteration  steps.  As  one 
can  see,  the  target  shape  and  localization  were  poorer 
them  those  in  the  high-contrast  case.  In  addition,  an 
artifact  appeared  at  the  edge  of  the  image. 

From  tile  coregistered  ultrasound  images,  we  ob¬ 
tained  spatial  distributions  of  the  two  targets  and 
specified  target  regions.  Figures  12(a)  and  12(c) 
show  the  —6-^  contour  plots  of  Figs.  li(a)  and  11(c). 
Applying  the  same  reconstruction  scheme  to  these 
specffic  regions,  we  obtained  Figs.  12(b)  and  12(c)  in 
one  iteration.  The  reconstruct^  absorption  coeffi¬ 
dents  were  0.2357  smd  0.219  cm”^  for  the  two  high- 
contrast  target  cases  emd  0.123  and  0.131  cm~^  for 
the  low-contrast  case.  We  can  see  much  better  im¬ 
provement  in  the  low-contrast  target  case  when  we 
compare  Fig.  12(d)  with  Fig.  11(d).  This  example 
demonstrates  that,  when  tiie  tarots  are  visible  in 
ultrasound  imeiges,  their  morphology  information 
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Elg.  si.  Reconstructed  NIR  images  for  shallow  targets  (LS  cm  in  depth,  1  cm  in  diameter,  and  the  fitted  background  |1b  and  m'  are  0.015 
and  5.36  cm~\  respectively).  The  left  column  corresponds  to  images  of  a  high-contrast  target  (pa  °  0.25  cm~*),  and  Uie  right  column 
corresponds  to  images  of  a  low-contrast  target  (p;,  •>  0.1  cm~^).  Each  row  is  related  to  a  specific  hole  size:  (a)  and  (b)  no  hole,  (c)  and 


(d)  2  cm  X  2  cm,  (e)  and  (f)  3  cm  x  3  cm. 


provided  by  ultrasound  can  be  used  to  guide  the  op¬ 
tical  reconstruction  in  the  specified  regions. 

The  result  regarding  the  iteration  step  is  signifi¬ 
cant.  As  we  discussed  above,  there  is  no  known 
stopping  criterion  to  terminate  the  iteration  because 
it  is  difficult  to  estimate  the  noise  level  in  the  mea¬ 
surements.  With  the  a  priori  target  depth  and  spa¬ 
tial  distribution  provided  by  coregistered  ultrasound, 
we  can  obtain  an  accurate  optical  absorption  coeffi¬ 
cient  in  one  iteration.  Therefore  no  stopping  crite¬ 


rion  is  needed  for  the  inversion  algorithms. 
However,  this  result  wiU  need  to  be  further  evaluated 
with  more  samples  of  different  contrasts. 

5.  Discussion 

Commercial  ultrasound  scanners  use  one¬ 
dimensional  probes  that  provide  2-D  images  of  x-z 
■views  of  the  tergets,  wherex  andz  are  the  spatial  and 
propagation  dimensions,  respectively.  Sudi  x—z  im¬ 
ages  cannot  coregister  with  NIR  images,  which  are 
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Table  2.  Parameters  o(  Reconstructed  images  for  StiaOow 
High-Contrast  (m,  *:  0.2S-cm"’)  and  Lotv-Contrast  (|i.  «  ai-cm”’) 
Targets* 


Probe  Type 

Parameter 

No  Hole 

2  cm  X  2  cm 

3  cm  X  3  cm 

High  contrast 

0.250 

0.194 

'  0.118 

{cm"‘) 

FWHM 

1.32 

1.61 

2.08 

(cm) 

Artifacts 

-7.98 

-12.7 

-9.76 

(dB) 

Low  contrast 

ita(max) 

0.100 

0.091 

0.042 

(cm"‘) 

FWHM 

1.88 

2.11 

3.17 

(cm) 

Artifacts 

-6.25 

-7.44 

-0.65 

(dB) 

*The  fitted  background  and  |l,'  are  0.015  and  S.3S  cm~\ 
respectively. 


obtained  from  x-y  views  of  the  taiigets.  Our  current 
2-D  ultrasound  array  is  capable  of  providing  x-y 
views  of  the  targets,  which  can  be  used  to  coregister 
with  NIR  images.  However,  the  array  is  sparse  and 
therefore  the  image  resolution  is  not  state  of  the  art. 
Nevertheless,  its  spatial  resolution  is  comparable  to 
NIR  imaging  and  can  be  used  to  guide  NIR  image 
reconstruction.  With  3-D  ultrasound  guidance,  only 
one  iteration  is  needed  to  obtain  accurate  absorption 
coefficients.  This  result  is  significant  because  no 
stopping  criterion  is  necessary.  More  studies  with  a 
variety  of  target  contrasts  and  locations  will  be  per¬ 
formed  to  verify  this  result. 

We  purchased  a  2-D  state-of-the-art  ultrasound  ar¬ 
ray  of  1280  transducer  elements  and  we  are  building 
a  multiplexing  unit  for  our  64-chaimeI  electronics. 
In  admtion,  the  new  2-D  transducer  size  is  approxi¬ 
mately  2  cm  X  3  cm,  which  is  in  the  neighborhood  of 
the  optimal  hole  size  we  foimd  through  this  study. 
With  the  new  2-D  ultrasound  transducer,  we  will  be 
able  to  obtain  hi^-resolution  ultrasound  images  and 
delineate  the  target  boundaries  with  finer  details  for 
optical  reconstruction. 

Ultrasound  contrast  depends  on  lesion  acoustic 
properties,  and  NIR  optical  contrast  is  related  to  le¬ 
sion  opticrd  properties.  Both  contrasts  exist  in  tu¬ 
mors,  but  the  sensitivities  of  these  two  modalities 
may  be  different.  It  is  possible  that  some  early- 
stage  cancers  have  NIR  contrast  but  are  not  detect¬ 
able  by  ultrasound.  It  would  be  desirable  if  we  could 
obtain  sensitivity  of  optical  imaging  alone.  How¬ 
ever,  light  scattering  is  a  main  problem  that  prevents 
the  accurate  and  reliable  localization  of  lesions.  It  is 
also  possible  that  some  lesions  have  acoustic  contrast 
but  no  NIR  contrast  or  low  NIR  contrast.  Currently, 
ultrasound  is  routinely  used  as  an  adjunct  tool  to 
x-ray  mammography;  the  combined  sensitivity  of 
these  two  modalities  in  breast  cancer  detection  is 


more  than  90%.^  Recently,  ultrasound  has  also 
been  advocated  to  screen  dense  breasts.®®  We  antic¬ 
ipate  that  our  combined  imaging  will  add  more  spec¬ 
ificity  to  the  ultrasonically  detected  lesions. 

In  the  reported  phantom  studies,  we  assigned  zero 
perturbations  to  &e  regions  where  no  targets  were 
present.  In  clinical  studies,  we  plan  to  segment  the 
ultrasound  images  and  spe<^  different  tissue  types 
as  weU  as  suspicious  regions  in  the  segmented  im¬ 
ages.  We  will  then  reduce  the  reconstructed  optical 
ufiJmowns  by  assigning  unknown  optical  properties 
to  different  tissue  types  as  well  as  to  suspicious  re¬ 
gions.  Finally,  we  reconstruct  the  reduced  sets 
of  unknown  optical  properties.  We  expect  a  more 
accurate  estimation  of  reconstructed  optical  proper¬ 
ties  and  fast  convergence  speed,  as  reported  in  this 
paper.  However,  it  is  still  too  early  to  judge  the 
clinical  performance  of  the  combined  method;  further 
clinical  studies  are  needed. 

Probing  regions  of  the  banana-shaped  diffusive 
photons  depend  on  source-detector  separations  and 
measurement  geometry.  For  a  semi-infinite  geome¬ 
try,  the  probing  regions  extend  further  into  the  me¬ 
dium  when  source-detector  separation  increases. 
This  is  why  we  have  multiple  source-detector  pairs  of 
various  ■  separations  to  detect  targets  at  variable 
depths  from  0.5  to  4  cm.  Of  course  it  is  hard  to 
adueve  uniform  sensitivity  in  the  entire  region  of 
interest.  For  example,  a  superficial  target  (~1  cm 
deep)  would  cause  strong  perturbations  when  it  is 
close  to  a  source  or  a  detector,  but  will  result  in  much 
weaker  signals  when  it  is  located  deeper.  Normal¬ 
ization  of  scattering  photon  density  waves  with  re¬ 
spect  to  the  incident  waves  makes  it  possible  for 
reconstruction  algorithms  to  handle  the  huge  dy¬ 
namic  range  of  signals  and  to  detect  a  target  as  deep 
as  4  cm.  This  normalization  procedure  was  applied 
to  the  reconstruction  algoritW  used  to  obtain  the 
reported  images. 

In  this  study,  the  target  absorption  coefficient  was 
reconstructed  from  the  measurements.  Because  the 
target  p.,'  was  similar  to  the  backgroimd  p.,',  the 
coupling  between  iig  and  p.,'  in  our  measurements 
was  negligible.  We  also  performed  experiments 
with  gel  phantom  made  witii  Intralipid  with  a  con¬ 
centration  similar  to  that  of  the  backhand  and  did 
not  observe  perturbation  beyond  the  noise  level. 
Similar  reconstruction  studies  can  be  performed  for 
scattering  coefficients  as  well.  Simultaneous  recon¬ 
struction  of  both  absorption  and  scattering  coeffi¬ 
cients  is  also  possible.  Because  the  eigenvalues  of 
the  absorption  and  scattering  weight  matrices  are 
significantly  different,  good  regulation  schemes  are 
needed  for  simultaneous  reconstruction.  This  sub¬ 
ject  is  one  of  our  topics  for  further  study. 

6.  Summary 

We  have  constructed  a  near-real-time  imager  that 
can  provide  coregistered  ultrasound  and  NIR  images 
simultaneously.  This  new  technique  is  designed  to 
improve  the  specificity  of  breast  cancer  dis^osis. 
Because  the  ultrasound  transducer  needs  to  occupy 
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Pig.  10.  Deep  target  (2.5  cm  in  depth,  1  cm  in  diameter)  of  low  optical  contrast  (|t„  =  0.10  cm"*  and  fitted  backgroimd  and  p.,'  are  0.02 
■n<l  6.08  cm"*,  respectively),  (a)  A-scan  line  of  the  reflected  ultrasound'  pulse-echo  signal  indicating  the  target  depth,  (b)  Absorption 
image  of  the  low-contrast  target  obtained  firom  optical-only  reconstruction,  (c)  Ultrasound-guided  reconstruction  at  target  depth. 


the  central  region  of  the  combined  probe,  a  series  of 
experiments  were  conducted  to  investigate  the  effects 
of  missing  optical  sensors  in  the  middle  of  the  com¬ 
bined  probe  on  the  NIR  image  quality-  Our  results 
have  shown  that,  as  long  as  the  central  ultrasound 
transducer  area  is  in  the  neighborhood  of  2  cm  X  2 


cm,  essentially  similar  reconstruction  results  as 
those  of  no  missing  optical  sensors  in  the  middle  of 
the  combined  probe  can  be  obtained.  This  2  cm  x  2 
cm  dimension  is  approximately  the  size  of  most  com¬ 
mercial  ultrasoimd  phased-array  transducers. 
When  the  central  missing  optical  sensor  area  is  in- 
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Elg.  11.  Simultaneously  obtained  liltrasound  and  NIR  absorption  images.  The  fitted  background  and  p,'  are  0.017  and  4.90  cm“  , 
respective^,  (a)  Ultrasound  and  (b)  NIR  absoiption  image  of  two  high-contrast  targets  (target  |i^  =  0.26  cm” ‘).  (c)  Ultrasound  and  (d) 
NIR  image  of  two  low-contrast  targets  (target  »=  0.10  cm~^).  In  both  high-  and  low-wntrast  cases,  the  two  targets  were  located  at  2.5 
cm  in  depth.  . 


creased  to  3  cm  X  3  cm,  however,  the  reconstructed 
values  are  obviously  lower  than  real  values.  If  we 
increase  the  iteration  steps,  artifacts  in  the  recon¬ 
structed  images  would  soon  become  dominant. 

With  the  target  3-D  distribution  provided  by  coreg¬ 
istered  ultrasound,  significant  improvements  in  algo¬ 
rithm  convergence  and  reconstruction  speed  were 
achieved.  In  general,  a  priori  target  depth  informa¬ 
tion  guides  the  inversion  algorithm  to  reconstruct  the 
heterogeneities  at  the  correct  spatial  locations  and 
improves  the  reconstruction  speed  by  an  order  of 
magnitude.  In  addition,  the  a  priori  target  spatial 
distribution  can  further  reduce  the  iteration  to  one 
step  and  also  obtain  accurate  optical  absorption  co¬ 
efficients.  Given  the  fact  that  no  known  stopping 
criterion  is  available,  this  result  is  significant  be¬ 
cause  no  iteration  is  needed.  However,  this  result 
will  need  to  be  evaluated  with  more  samples  of  dif¬ 
ferent  contrasts. 
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